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membrane lipids and in the double layers lead to corrections to the membrane elastic 
moduli which are analyzed here. We show that, especially in the low salt limit, i) 
the electrostatic contribution to the membrane's surface tension due to the Debye 
layers crosses over from a quadratic behavior in the externally applied voltage to a 
linear voltage regime, ii) the contribution to the membrane's bending modulus due to 
the Debye layers saturates for high voltages. Nevertheless, the membrane undulation 
instability due to an effectively negative surface tension as predicted by linear Debye- 
Hiickel theory is shown to persist in the nonlinear, high voltage regime. 
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1. Introduction 

Bilayers formed from lipid molecules are an essential component of the membranes 
of biological cells. The mechanical properties of membranes at equilibrium are 
characterized by two elastic moduli, the surface tension and the bending modulus pQ. 
These moduli typically depend on electrostatic properties, and their modifications in 
case of charged membranes/surfaces in an electrolyte have been examined theoretically 
in the period 1980-90's as reviewed e.g. in Ref. [2]: they have been first derived 
by Winterhalter and Helfrich [3] within linearized Debye-Hiickel (DH) approximation, 
then by Lekkerkerker jl] in the nonlinear Poisson-Boltzmann (PB) regime for charged 
monolayers and by Ninham et al. |5j for charged symmetric bilayers. Later on Helfrich 
et al. revisited the question of the electrostatic corrections to the bending modulus of 
charged symmetric bilayers [6]. 

Nowadays, the study of deformations of membranes or vesicles in electric fields is 
an active field of research linked to many biotechnological applications. For instance, 
the application of electric fields is used to produce artificial vesicles from lipid films 
(electroformation), or to create pores in vesicles (electroporation), which is an important 
route for drug delivery. Both processes are widely used experimentally although they 
are not well understood theoretically. The effects of electric fields on giant unilamellar 
vesicles have been reviewed recently in Ref. [Tj. This system shows a rich panel of 
possible behaviors and morphological transitions depending on experimental conditions 
- electric field frequency, conductivities of the medium and of the membrane, salt 
concentration, etc. These observations are supported to a large extend by theoretical 
modeling [8]. 

Originally it was Helfrich et al. who pointed out that the deformation of lipid 
vesicles in electric fields can serve as a means to determine the electrostatic corrections 
to the membrane elastic moduli [9]. Besides this observation on the macroscopic scale 
other techniques can provide valuable insights into the moduli corrections such as 
AFM, impedance spectroscopy [ID] , neutron reflectivity pZE| and X-ray scattering [12] , 
Recently [12], X-ray scattering experiments have been carried out on a system of two 
superposed lipid membranes in an AC electric field. In trying to analyze this data, we 
noted that these experiments have been carried out at relatively high voltages, in a 
regime where the linearized DH approach may not be applicable. In order to describe 
such a situation theoretically, we extend previous work [T3J [HI [15] based on the DH 
approach to the nonlinear PB regime, which is more suitable for realistic situations in 
which the induced surface charges on the membrane are large. 

In this paper, we present a simple approach to calculate electrostatic corrections in 
the elastic moduli of a quasi-planar lipid membrane. The membrane is assumed to be 
non-conductive to ions, non-permeable to water and electrically neutral, it is subjected 
to a normal DC electric field and embedded in an electrolyte described by the Poisson- 
Boltzmann equation. In this situation, the electric field leads to an accumulation of 
charges on both sides of the membrane, which affect the mechanical properties of the 
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membrane. The electrostatic corrections to the elastic moduli can be used to analyze the 
instability of a lipid membrane in an applied DC electric field. In contrast to previous 
work, based on a free energy approach [H El [16] , our method is based on a calculation 
of the stress (or force) balance at the membrane surface using electrokinetic equations 
[T7] . This method is related to the work of Kumaran [18] who used a similar approach 
in the context of equilibrium charged membranes. Two points are worth mentioning: 
first, our approach is able to describe the capacitive effects of the membrane and of the 
Debye layers while keeping the simplicity of the zero thickness approximation on which 
most of the literature on lipid membranes is based. Second, our approach can include 
non-equilibrium effects which can not be described within the free energy approach. For 
instance, in Refs. [131 El US] we investigated the effects of ionic currents flowing through 
the membrane, which in turn affect the fluid flow near the membrane. Other types of 
non-equilibrium effects that could be included in that framework are those arising from 
the stochasticity of ion channels. 



2. Model equations 

We consider a steady (DC) voltage V between two electrodes at a fixed distance 

z = ±L/2, applied to an initially flat membrane located at z — 0. The membrane is 

embedded in an electrolyte of monovalent ions with densities n + and n~ . The membrane 

is treated as non- conductive for both ion species and is (effectively) uncharged; thus we 

focus solely on capacitive effects. A point on the membrane is characterized within 

the Monge representation by a height function h(r±), where r_i_ is a two-dimensional 

in-plane vector in the membrane. 

In the electrolyte, the electric potential <fi obeys Poisson's equation 

1 2 
V 2 </> = — (en + — en~) = — p. (1) 

Here e is the elementary charge, e is the dielectric constant of the electrolyte and we 
have introduced half of the charge density, 

n + — rT 

P = e • (2) 

For the sake of simplicity, we assume a symmetric 1 : 1 electrolyte, so that far away 
from the membrane n + = rT = n*, and the total system is electrically neutral. 
The densities of the ion species obey the Poisson-Nernst-Planck equations 

«9 t n ± + V-j ± = 0, (3) 

with ionic current densities 

± _ n / w„± -r- „± e 



3-D I -Vn* t n ^V0j , (4) 

where ksT is the thermal energy. For simplicity we consider the case where both 
ion types have the same diffusion coefficient D and neglected various corrections for 
concentrated solutions |19|. 
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As boundary conditions (BC), the potential at the electrodes is externally held at 

-4) -4- 

This BC is oversimplified for real electrodes, but captures the main effects of the electric 
field, see the discussion in Ref. [T5]. The distance between the electrodes is assumed to 
be much larger than the Debye length, L ^> \r> = k,^ 1 . In that case, the bulk electrolyte 
is quasi- neutral, n + = n~ = n*, with negligible charge density (compared to the total 
salt concentration), so that far from the membrane 

p(z = ±%)=0. (6) 



2, 

The BC at the membrane is crucial to recover the correct physical behavior. Let n 
be the unit vector normal to the membrane. We use the Robin-type BC 

A m (n • V)(j)\ z=h + = A m (n • V)<f)\ z=h - = (p(h + ) - <f>[h~) , (7) 

where 

A m = — d (8) 

Cm 

is a length scale that contains the membrane thickness d and the ratio of the dielectric 
constant of the electrolyte, e, and of the membrane, e m . This BC was originally 
developed for electrodes sustaining Faradaic current [201 121] or charging capacitively 
[22] . In Refs. [231 [HI [15] , this BC was derived, and was shown to properly account for 
the jump in the charge distribution which occurs near the membrane as a result of the 
dielectric mismatch between the membrane and the surrounding electrolyte. 

3. Poisson-Boltzmann approach for a membrane in an external potential 

Here we show how the well-known solution of the Poisson-Boltzmann (PB) equation for 
a single charged plate in an electrolyte can be used to describe the present situation of 
a capacitive membrane with induced Debye layers in an external potential. 

In a steady state situation and when there is no electric current through the 
membrane, one obtains from Eqs. ([3j H} 

- Vn ± ^n ± -^-V0 = O. (9) 
After a direct integration using the BCs from above, one obtains 

n± = n*e T ^ (2) "^ , (10) 
and insertion in Poisson's equation then yields the Poisson-Boltzmann (PB) equation 

Linearization (for — 1) leads to the well-known Debye-Hiickel (DH) equation, 

V 2 = L(z) - ^\ , (12) 
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where 

, 2e 2 n* , . 

and K~ x = is the Debye length that defines the characteristic length scale for charge 
relaxation in the electrolyte. 

The nonlinear PB equation (fTTI) for the planar case can be integrated analytically 
[2], leading to 



e \ 1 — ce~ KZ j 2 



-KZ \ 2 



) , (15) 



1 =F ce 

for z > 0. The expressions for z < can be obtained using the symmetry of the system: 
4>(—z) = —<f>(z), n (— z) = —n (z). 

Just as in the classical PB solution for a single charged plate in contact with an 
electrolyte, the non-dimensional parameter c is determined by the BC for the potential 
at the membrane. For a flat charged surface surrounded by an electrolyte [2], c is given 
by a simple quadratic equation and can be expressed in terms of the ratio of the two 
characteristic length scales: the Gouy-Chapman length b = 2e e ^ T (with |er| the charge 
density of the surface) and the Debye length Ad = k -1 - In contrast, in the case of a 
membrane in an electric field, we obtain from Eq. (J7|) the following nonlinear equation 

4KA m ^ + -^-=41n(l±^) . (16) 



c 2 -l k B T \l-c, 
Note that two dimensionless parameters enter this equation: the ratio of electrostatic 
to thermal energy, 

and the dimensionless parameter [21] 

A m — KA m — — — — — — — , [LtS) 

which quantifies the electrical coupling between the membrane and the Debye layers. 
More precisely, depending on whether this parameter is large or small with respect to 
one, the capacitance of the diffuse part of the double layer, Cp = en, or the capacitance 
of the membrane, C m = e m /d, dominates the overall voltage drop. 

The non-dimensional parameter c = c(A m , V) which is determined by Eq. ([IS]) 
is related to the potential at the membrane, 4>{0 + ), and to the charge density at the 
membrane, p(0 + ), by the following relations 

0(0+) = ^-artanh(c) + — , (19) 

p(0+) = en* sinh [4 artanh(c)] . (20) 

From Eq. (fl9|) . one can see that the values of c must be restricted to the interval [0, 1]. 
In the limit of small voltages, V" <C 1, there is a linear relation between c and the charge 
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Figure 1. a) Solution of Eq. (fT6|) for different amounts of salt: k = 10 6 m _1 (pure 
water; black), k = 10 7 m -1 (red), k = 10 8 m _1 (green), k = 10 9 m _1 (blue), b) Charge 
distribution (in units of Cm -3 ) for z > : comparison of the nonlinear solution (blue) 
with the linearized Debye-Hiickel solution (green). The asymmetric distribution of 
positive (red dashed) and negative charges (black dashed) for the nonlinear solution is 
also shown. Parameters used: e = 80eo (water) and e m = 2eq for the membrane. As 
the membrane thickness is typically d = 5nm this leads to A m = j-d = 200nm. For 
b) we used k = 10 6 m _1 , corresponding to en* ~ 9.16 Cm~ 3 , and V — 0.1V. 



distribution at the membrane, since c = = z — i^rr-^r, in accordance with the 

4en* 8era*(2+KA m ) ' 

calculations based on the DH approximation of Ref. [15] . 



Fig. [Tt) shows the solution c of Eq. ( fT6l) as a function of external voltage for various 
values of n, which corresponds to varying the amount of salt since n oc \/n*. Clearly, 
for low salt, the linear approximation remains valid only for rather small voltages (for 
instance it holds only for V < 0.1V in pure water), while for high salt V = 5V is 
still in the linear regime. Fig. [Tb) shows a comparison of the charge distribution for 
the nonlinear PB (blue) and the linear DH solution (green). As expected, the figure 
shows that the DH approximation underestimates the surface charge on the membrane 
layers as compared to the PB calculation. The figure also shows the distribution of the 
positive and negative ions, which both tend to n* far from the interface as a result of 
electroneutrality. 

Although the present situation differs from the case of a single charged plate in 
an PB electrolyte, the structure of the solutions Eqs. ( IT4] [T5|) is very similar in both 
problems. Because of this, there is an equivalent to the Contact theorem [2], which 
relates in the single charged plate problem the surface charge density to the limiting 
value of the potential/ionic density at the plate: namely, one can give the effective 
surface charge a c s for the charged plate problem that creates the same volt age/ charge 
distributions as the capacitive membrane in the external field. This effective surface 
charge reads 



of ff = Aen*k B T 



cosh ( ^(0(0+) -V/2) 



(21) 
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4. Corrections to membrane elastic moduli 

Surface tension. The electrostatic corrections to the surface tension can be calculated 
directly from the stresses acting on the membrane in the flat configuration (also called 
the base state) as explained in Refs. [24] US]- The total stress tensor reads 

r l3 = -pSij + v (diVj + d^) + e (E i E j - ^<%# 2 J (22) 

which contains the pressure, the viscous stresses in the surrounding fluid (the electrolyte) 
and the Maxwell stress due to the electrostatic field, rj is the viscosity of the electrolyte 
and v its velocity field. The electric field is given by E = — V0. 

In the base state, the electric field is oriented in ^-direction, and the condition 
V • t = implies that d z p = |<9 Z ((<9 Z 0) 2 ) = —2pd z <f). Using Eqs. (]T4"] [15]) and imposing 
p{z — > oo) = 0, this is readily solved by 

p 2ft2 

p(z > 0) = lQn*k B Tc 2 g ? , (23) 

(c + e KZ ) (c — e KZ ) 

and similarly with z — > — z for z < 0. 

Let us call 5 a closed surface englobing the membrane with the normal vector field 

n. The force acting on the surface S in the x-direction (chosen to be the direction of 

the lateral stress) can be calculated from the stress tensor as F x = J s x • r ■ n dS. Since 

r is divergence free, the surface S can be deformed, for convenience to a cube of size L, 

and it is easy to see that the integral is non-zero only on the faces of the cube with the 

normal along ±x. With dS = Ldz and AS = F x /L for n = +x, we arrive at 

nL/2 

AS= / r xx {z)dz. (24) 

J-L/2 

From Eq. (122]) . one has t xx {z) = —p(z) — |(<9 2 0) 2 , where <j>(z) and p(z) are the potential 
and the pressure in the base state given above. Upon integration (using Lk ^> 1), one 
obtains the corrections to the surface tension, as a sum of two terms. 
First, there is the external contribution due to the Debye layers, 
32n*k B T c 2 

K 1 — C 2 

Second, there is the internal contribution due to the electric field inside the membrane 
(cf. Ref. pj)]), which is given by E™ = — \ (4>(0 + ) — 0(O~)). That correction to the 
surface tension is AE m = — e m c?(i?™) 2 , or explicitly 

AE^-^LiM^fl+fW] 2 . (26) 
a e \1 — c J 

Note that both corrections to the surface tension are negative, which means that 

these corrections can lead to an instability as soon as the total surface tension S to4 = 

So + AS K + AS m (which is the sum of the bare tension Eo plus the above corrections) 

becomes negative [2"B] . 

Bending modulus. To obtain the correction to the bending modulus we perform, 

as detailed in Ref. [15J for the DH calculation of the potential at first order in 



AE K = _ B _____ (25) 
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the membrane height. Then, by solving the hydrodynamics problem of the electrolyte 
around the membrane (in Stokes approximation), one determines the pressure and 
obtains the total stress tensor. The growth rate of membrane fluctuations, s{k±), where 
\t± is the wave vector in the membrane plane, is then determined by imposing that the 
discontinuity of the normal-normal component of the total stress tensor at the membrane 
has to equal the membrane restoring force: 



- [(r zz ,i(z = 0+)-t ZZj1 (z = 0-)] 
Here Fh is the standard Helfrich free energy, 



dh(r ± ) ' 



H 



(27) 



(2f 



and S and K are the bare surface tension and the bare bending modulus of the 
membrane, respectively. Expanding the left hand side of Eq. ( J27|) in powers of yields 
the growth rate of the form 



rjk±s(k_i 



AS. 



\(K 



AK K + AK m ) k\ 



(29) 



Details of the calculations can be found in Appendix A The surface tension 



corrections calculated above in Eqs. (125] 126]) can be recovered by this method, which 
provides a self-consistency check. For the bending modulus one obtains 



AK K 



AK n 



8n*k B Tc 2 (3 



K 



1 + C 2 



_d*_ 
12 + Elen 



2k B T c(l - c 2 ) 



d 



(30) 
(31) 



1 + c 2 

for the external contribution due to the Debye layers and the internal contribution due 
to the voltage drop at the membrane, respectively. The field inside the membrane is 
given by E™ = - J [-^f- In (±±^) + V] . 



5. Discussion 



Let us now discuss the nonlinear electrostatic effects on the membrane elastic moduli in 
the limits of low and high applied voltages. 

Low voltage regime. In the low voltage limit, F< 1, a solution of Eq. ffl6|) to 
linear order in c, yields 

4(2 + A m ) A(2 + \ m K)k B T 4en* ' v ; 

Here p(0 + ) is half the charge density at the membrane, corresponding to the quantity 
called p m in Ref. jl5j for the DH case. Expanding <j)(z) and p(z), as well as the corrections 
to the moduli AS K , AE m , AK K and AK m for small c, one exactly recovers all of the 
results given in Ref. [15]. Specifically, all corrections to the moduli scale quadratically 
with the external voltage, ex V 2 . This is due to the fact that both the potential and the 
induced charge are proportional to the applied voltage. 
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Figure 2. Corrections to the surface tension for k = 10 m (pure water) as a 
function of voltage. Corrections from the Debye layer are in red, corrections due to 
the voltage drop at the membrane are in blue, a) Nonlinear PB result. One sees that 
the correction due to the Debye layer scales linearly in the voltage for high voltages 
(cf. the black dashed straight line), while the contribution from inside the membrane 
scales quadratically in the voltage, b) Comparison of the nonlinear PB result (solid 
lines) with the linear DH result (dashed lines). 



High voltage regime. In the opposite limit, V ^> 1 implies c — > 1. Introducing 
a = 1 — c and expanding Eq. ( IT6l) for small a, one gets — ^ L +Alna + V+X m — 4 In 2 = 0. 
For high values of V one can neglect the last two (constant) terms and gets 

(33) 



W (±A m e^) V/4 ' 

where W(x) is Lambert's function, i.e. y = W(x) is the solution of ye y = x. As a 
result, the nonlinear electrostatics strongly effects the corrections to the moduli from 
the Debye layers. In fact, the external contribution to the surface tension scales as 

AE ^ >>1)K __^_Z_ (34) 

instead of A£ K (V <C 1) oc — , , 2 : in the high voltage regime, the external surface 
tension correction scales linearly with the voltage instead of quadratically. The crossover 
from V 2 to V is salt dependent, i.e. depends on the value of A m = A m K, cf. Fig. [TJ In 
contrast, the external surface tension correction AE m remains a quadratic function of 
the applied voltage. 

Fig. [2] displays both corrections to the surface tension as a function of voltage for 
low salt (pure water). Fig. Eh) shows that due to the crossover to a linear voltage 
dependence for the external surface tension correction, rapidly AE m (the blue curve) 
wins and AE K (red curve) becomes negligible for V > 1.5 V. Fig. Eb) shows a comparison 
of the nonlinear PB result (solid lines) with the linear DH result (dashed lines). For 
V < 0.2V, there is agreement and the external contribution dominates. However, for 
higher voltages the linearized DH solution becomes completely misleading: it predicts 




Figure 3. Corrections to the bending modulus for k = 10 6 m _1 (pure water) as 
a function of voltage, a) Corrections from the Debye layer, b) Corrections due to 
voltage drop at the membrane. Nonlinear PB calculations are solid, linearized DH 
calculations are dashed (note the different scales on the left and right panels). 



that |A£ K | > |AE m | for all voltages (and both proportional to V 2 ), while already above 
V > 0.4V" the internal contribution exceeds the external contribution. 

For the bending modulus corrections, the differences between the DH and PB 
models are even more striking: As a function of the applied voltage, the external 
contribution levels off at a constant value 

» 1) = ^L^l _ ™ . (35 , 

K a 1 + C z K 

In contrast, the internal contribution AK m continues to grow quadratically in the high 
voltage limit. Fig. [3^) and b) show comparisons of the nonlinear PB solutions (solid) 
and the linear DH solutions (dashed) for both contributions to the bending modulus. 
Note, however, the different scales: Although the external contribution saturates, see 
Fig. [HK), this value is larger by two orders of magnitude than the still growing value 
from the internal contribution, cf. Fig. [3b). In conclusion, for pure water and voltages 
~ 0.2 — 5V, the total correction to the bending modulus will appear constant within this 
voltage interval. Only for still higher voltages the quadratic growth due to the internal 
contribution AK m will dominate. 

Membrane instability. Let us briefly discuss the consequences for the membrane 
undulation instability. As mentioned above, an instability develops in this system when 
the total surface tension, the sum of the bare value and the electrostatic corrections, 
becomes negative. The threshold value for the voltage, V c , has been calculated in 
Ref. [15] within linearized DH theory and reads (in case of a non- conductive membrane) 

y2 = Sprf (2 + K\ m ) 

k (ne m \ 2 m + ed) ' 

This curve is shown in Fig. H] as the blue line. The red line, in contrast, shows the 
nonlinear PB result given by numerical solution of = S tot = S + A£ K + AX m 
with Eqs. (121)1 for the surface tension corrections and c(A m , V) calculated from 




Figure 4. Threshold voltage V c for the membrane undulation instability as a function 
of salt (i.e. k). Nonlinear PB result (red curve) vs. linear DH result (blue curve). 
Parameters: k — 10 6 m _1 (pure water); bare surface tension So = 10 _3 Nm _ ; 
X m = 200nm. 



Eq. (fT6|) . In the high salt limit, the Debye layers shrink to zero and as a result the 
external contribution vanishes. For the internal contribution, the inside field is exactly 
calculated from the DH approach, as can be seen from the behavior of the parameter c 
in Figure dt. Thus, both curves merge and attain the limiting value given by Eq. (I36p . 
VJk — > oo) = \ In contrast, for low salt, the linear result overestimates the 
threshold value, since it underestimates the induced charges at the membrane, as shown 
in Fig. [lb). 

In addition to the threshold voltage for the instability, the most unstable wave 
number associated to the instability will be affected by the electrostatic nonlinearities 
of the PB equation. This wave vector corresponds to the maximum of the growth rate 
given by Eq. ( 1291) with respect to k±. In the low salt regime, since the system is more 
unstable in the presence of electrostatic nonlinearities, cf. Fig. HI in general the wave 
vector will be increased as compared to the predictions from the DH approach. 

6. Conclusions 

We have investigated the nonlinear electrostatic effects of an external DC electric field 
on a purely capacitive membrane, which is non-conductive for the ions and bears no 
fixed charges, in an electrolyte. We have calculated in the nonlinear Poisson-Boltzmann 
regime the corrections to the membrane elastic moduli - both the external ones due 
to the Debye layers surrounding the membrane and the internal corrections due to the 
electric field inside the membrane. Strong deviations from the linear Debye-Hiickel 
behavior have been found in the low salt regime at already rather moderate voltages. In 
particular we have shown that the external contribution to the surface tension crosses 
over from a quadratic dependence on the externally applied voltage as predicted by the 
linearized theory to a linear voltage dependence. In contrast, the internal contribution 
remains quadratic and becomes dominating at high voltages. The external contribution 
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to the bending modulus even saturates for high voltages, while the internal contribution 
remains quadratic in voltage. 

In addition, our work confirms that the surface tension still grows in absolute value 
with the voltage, which means that the membrane undulation instability present in 
the DH theory (due to an effectively negative surface tension) persists in the nonlinear 
PB regime. The nonlinearities affect the threshold in voltage and the characteristic 
wavelength of the instability. 

The method presented here can serve as a starting point for further extensions. 
One example would be to include fixed charges in addition to induced charges, similarly 
as studied in [26]. Another possible direction would be to improve the description of the 
membrane, cf. the recent work in Ref. [27J where thickness fluctuations of the membrane 
and fluctuations of the lipid dipole orientations within the membrane are accounted for in 
a comprehensive continuum model for a membrane in a normal DC electric field. Other 
possible extensions could include other types of non-linear effects, for example due to 
membrane elasticity, due to inclusions of proteins such as ion channels or pumps in the 
membrane [28] , and also various relevant non-equilibrium effects, coupling electrostatics 
and hydrodynamics as in induced charge electro-osmosis [29]. Finally, in the biological 
context, the heterogeneity of the bilayer composition is another feature which is beyond 
the present model and which is likely to be important. 
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Appendix A. Details of calculation of bending modulus 

As this is the extension to nonlinear electrostatics of earlier work |15j , we only sketch the 
calculations. To solve the electrostatics problem to first order in the membrane height, 
one linearizes in h by writing 

0(k ± , z) = <f) {z) + 0i(k ± , z) , 
p(k ± , z) = p (z) + pi(k ± , z) , 
ra^ki, z) = tiq(z) + nf (kj_, z) , 

where (f>o, n$ are the base state solutions (flat membrane) as given in the main part and 
Po — f ( n o ~~ n o )• Quantities with subscript 1 are of order h. We used the definition of 
the Fourier transform for the in plane vector rj_, f(k±,z) = J dr±e~ l]iJ -' rj - f(r±, z). 

Assuming zero current through the membrane, the Poisson-Nernst-Planck equation 
linearized in h has solutions n\ = —n^i and n{ = n^ipi. Insertion into the PB equation 
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1 + 6c 2 e~ 2KZ + c 4 e" 4 ^) 



yields to linear order in h 

(Si - hi) 0! = K 2 

V 2 ±J (1 - c 2 e" 2 « 

which is solved (with BC <pi(z —> oo) = 0) by 



(A.l) 



.4 



\ 



4 



K 



^>2 gS/^ z 



k 2 ( e 2nz _ c 



2\ 2 



p 2kz _ ( 1 -k) 2 r .2 
e — c 



-Iz 



(A.2) 



where we have introduced I 2 = k\ + k 2 (note that for c « 1 one regains a simple 
exponential, 0i = A'e~ lz , as in Ref. [ISj). The constant A can be obtained from the 
BC, Eq. (CD, at order h 

^[(d^oj^-h+id^)^] =0 1 (O + )-0 1 (O-). (A.3) 

As by symmetry (f>(z) = —<p(—z) and 0i oc /i, one has 4>i(z) = (f)x(—z) and Eq. ( 1A.3I) 
simplifies to (^0i)| z=o = — (<9 2 0o)| z= o ' ^- ^ is thus independent of A m . The full 
expression for A is not needed, since later on we expand in k±. 

Since we study the case without ionic current through the membrane, the 
hydrodynamics problem (cf. Ref. [15]) around the membrane is trivial and one gets 

v z = h{k ± )s{k ± ) (1 + k ± z) e~ k ^ z (A.4) 

for the normal component of the velocity, as induced by a pure membrane bending mode 
[30] . Here s(k±) is the growth rate of the membrane fluctuations. The pressure is given 
by (in incompressible Stokes approximation) 



P 



-rjd z v z + 



k f 



(A.5) 



Ah.2 1 l2 
v rv ftj | 

Herein enters the bulk force due to the electric field acting on the charge distribution, 
which reads f = — 2pV0 = — 2p o V0i — 2piV(pQ + 0{h 2 ) with perpendicular component 

fj_ = - 2po(z)zkj>i(ki, z) . (A.6) 

The total normal stress at the membrane to linear order in h is, cf. Eq. (J2j 



T zz x 



-p + 2 V d z v z + -(d z( j ) j z -^-(d z 



(A.7) 



We now have to calculate the normal stress discontinuity at the membrane up to fourth 
order in k± and to balance it with the membrane restoring force. With the abbreviation 
[/]o = /(0 + ) — /(0~) the stress discontinuity can be written as 

[r zz ,i] z=0 = - [p] z=0 + 2 V [d z v z } z=0 + e [(&0o)(&0i)L=o - e m [(d z ^)(d z <PT)} z=±d/2 . (A.8) 
The balance reads 

dF H 



[t zz ,i] z =o 



dh(r ± ) 



and yields the growth rate of membrane fluctuations s(k±) of the form 



nk±s{k ± ) = ~ (E + AE) k\ 



{K + AK) k 



± ■ 



(A.9) 



(A.10) 



A Poisson-Boltzmann approach for a lipid membrane in an electric field 



14 



A little care has to be taken for the correction due to the internal field, cf. Ref. 
Note that we accounted for the finite membrane thickness d in Eq. flA.8j) above. As the 
internal field is constant, and due to symmetry, one can write 

e m [{d z ^){d z 4>T)]\ z=+d ,2 = ~2e m £o, m {d z 4>?\ z=+d/2 . (A.ll) 

The potential inside the membrane can be written (using again the symmetry, as well 
as p = inside; cf. also Ref. [H] for details) 

<t>?(k±, z) = <PT g ) ^T^i ( ek±Z + ^) ■ (A.12) 

4>T \k±j |) can be calculated approximately by using the outside solution, Eq. (IA.2j) . 
and imposing the BC at the membrane, leading to 

l, fj = <Pi (V, ~ h(k ± ) {d,ffl - d z <p ) lz=d/2 . (A.13) 

On the right hand side, now all quantities are known. Finally, one expands Eq. (1A.11I) 
and all the other quantities entering the total stress discontinuity, Eq. (1A.7[) . up to 
fourth order in k±. From Eq. (I A . 9 [) one then obtains the growth rate s(k±) of membrane 
fluctuations. 



References 



[1] U. Seifert, Adv. Phys. 46, 13 (1997). 

[2] D. Andelman in Handbook of Biological Physics, edited by R. Lipowsky and E. Sackmann, Elsevier, 

Amsterdam (1995). 
[3] M. Wintcrhalter and W. Helfrich, J. Phys. Chem. 92, 6865 (1988). 
[4] H. N. W. Lekkerkerker, Physica A 159, 319 (1989). 
[5] D. J. Mitchell and B. W. Ninham, Langmuir 5, 1121 (1989). 
[6] M. Winterhalter and W. Helfrich, J. Phys. Chem. 96, 327 (1992). 

[7] R. Dimova, N. Bezlyepkina, M. Domange Jordo, R. L. Knorr, K. A. Riske, M. Staykova, P. M. 

Vlahovska, T. Yamamoto, P. Yang and R. Lipowsky, Soft Matt. 5, 3201 (2009). 
[8] P. M. Vlahovska, R. S. Gracia, S. Aranda-Espinoza and R. Dimova, Biophys. J. 96, 4789 (2009). 
[9] M. Kummrow and W. Helfrich, Phys. Rev. A 44, 8356 (1991). 
[10] G. Wiegand, N. Arribas-Layton, H. Hillebrandt, E. Sackmann and P. Wagner, J. Phys. Chem. B 

106, 4245 (2002). 

[11] S. Lecuyer, G. Fragneto and T. Charitat, Eur. Phys. J. E 21, 153 (2006). 
[12] Linda Malaquin, These, Universite de Strasbourg (2009). 

[13] D. Lacoste, M. Cosentino Lagomarsino and J. F. Joanny, EPL 77, 18006 (2007). 

[14] D. Lacoste, G. I. Menon, M. Z. Bazant and J. F. Joanny, Eur. Phys. J. E 28, 243 (2009). 

[15] F. Ziebert, M. Z. Bazant and D. Lacoste, Phys. Rev. E 81, 031912 (2010). 

[16] T. Ambjornsson, M. A. Lomholt and P. L. Hansen, Phys. Rev. E 75, 051916 (2007). 

[17] A. Ajdari, Phys. Rev. Lett. 75, 755 (1995). 

[18] V. Kumaran, Phys. Rev. E 64, 011911 (2001); R. Thaokar and V. Kumaran, Phys. Rev. E, 66, 

051913 (2002); V. Kumaran, Phys. Rev. Lett., 85, 4996 (2000). 
[19] M. Z. Bazant, M. S. Kilic, B. D. Storey and A. Ajdari, Adv. Colloid Interface Sci. 152, 48 (2009). 
[20] E. M. Itskovich, A. A. Kornyshev and M. A. Vorotyntsev, Phys. Status Solidi A 39, 229 (1977). 
[21] M. Z. Bazant, K. T. Chu, and B. J. Bayly, SIAM J. Appl. Math. 65, 1463 (2005); K. T. Chu and 

M. Z. Bazant, SIAM J. Appl. Math. 65, 1485 (2005). 
[22] M. Z. Bazant, K. Thornton and A. Ajdari, Phys. Rev. E 70, 021506 (2004). 



A Poisson-Boltzmann approach for a lipid membrane in an electric field 



15 



[23] M. Leonetti, E. Dubois- Violette and F. Horrible, Proc. Natl. Acad. Sci. U.S.A. 101, 10243 (2004). 
[24] J. S. Rowlinson and B. Widom, Molecular Theory of Capillarity, Oxford University Press, Oxford 
(1982). 

[25] P. Sens and H. Isambcrt, Phys. Rev. Lett. 88, 128102 (2002). 

[26] P-C. Zhang, A. M. Keleshian and F. Sachs, Nature 413, 428 (2001). 

[27] R. J. Bingham, P. D. Olmsted and S. W. Smye, submitted to Phys. Rev. E (2009). 

[28] M. D. El Alaoui Faris, D. Lacoste, J. Pecreaux, J.-F. Joanny, J. Prost and P. Bassereau, Phys. 

Rev. Lett. 102, 038102 (2009). 
[29] M. Z. Bazant and T. M. Squires, Phys. Rev. Lett. 92, 066101, (2004). 
[30] F. Brochard and J. F. Lennon, J. Phys. (Paris) 36, 1035 (1975). 



